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Abstract 

The light-scattering spectrum for a Maxwell fluid in a steady state due to the 
presence of an external temperature gradient is computed. In such a fluid heat 
conduction is assumed to be governed by the classical Fourier law. The calculation 
is carried out through the use of fluctuating hydrodynamics. The effect of the 
non-Newtonian character of the fluid in the resulting spectrum is discussed. 
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1 INTRODUCTION. 

When light passes through a fluid, correlations over space and time in the density fluc- 
tuations are responsible for the scattering of light. In fact, the dynamic structure fac- 
tor S(\t,uj), which is the space-time Fourier transform of the space-time correlations of 
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density fluctuations in the fluid [Q], reflects the strength of the exchanges of energy and 
momentum between a light beam and the fluid as a function of wavevector k and fre- 
quency uj. For a one-component Newtonian fluid in equilibrium, the theory for S(k,L)) 
was first given correctly in 1934 by Landau and Placzek|| and later rederived using 
fluctuating hydrodynamics H and the approach advocated by Mountain [P],@,P|. In 
fluctuating hydrodynamics, random fluxes are added to the hydrodynamic equations (lin- 
earized about equilibrium) which then appear as stochastic equations of the Langevin 
type. The fluctuation-dissipation theorem provides a connection between the spectrum 
of spontaneous fluctuations in thermal equilibrium and the transport properties of the 
system. In contrast, the dynamic structure factor in Mountain's work is derived directly 
from the linearized hydrodynamic equations without resorting to any kind of fluctuating 
term added to the constitutive equations. It is worthwhile to note here that both routes 
lead to the same results and that these results have been confirmed many times by light 
scattering measurements in various systems. 

For nonequilibrium fluctuations in Newtonian fluids, the theory was developed rather 
more recently]?]] - |1| and the few available experiments seem also to confirm its validity [fT~9f| - 



1 21]. The extension of fluctuating hydrodynamics to steady states is made on the basis 
of the simple hypothesis that the fast degrees of freedom behave as "if they were in local 
equilibrium at the value of the thermodynamic parameters determined from the solution 
to the steady state problem" [fuH . With this statement it is assumed that the statistics 
of the stochastic forces behaves as if they were within a local equilibrium environment 
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around each point r of the system \T%\. It is usual under this hypothesis to rewrite the cor- 
relation functions of the stochastic fluxes by substituting the equilibrium thermodynamic 
properties and transport coefficients by their space dependent values. Another way to see 
this assumption is based on the Markov property: the random fluxes are not affected by 
the applied gradients since the correlation length is of microscopic order ] |T^|j . 

Studies of fluctuations in non-Newtonian fluids are also relatively abundant. Some of 
these systems have also been described under the scope of linear Langevin approaches, 



as is the case of viscoelastic fluids near equilibrium states |22| , suspensions in viscoelastic 
fluids under density and velocity gradients [I7| and colloidal suspensions including memory 



effects 18]. On the other hand, there are many treatments a la Mountain in this kind 



of fluids as well |23| - f27fl . Among these, although at the time meant in a different 
context, we want to mention here explicitly the study of (equilibrium) Rayleigh-Brillouin 
scattering in polyatomic fluids within an extended irreversible thermodynamic treatment 



1 24j carried out by Lopez de Haro, Rodriguez and Garcia-Colin a few years ago [25|. There, 
no explicit emphasis on the non-Newtonian character of the fluid was made, although of 
course the evolution equations for the trace and the symmetric traceless part of the stress 
tensor clearly imply this character. More recently, Vazquez and Lopez de Haro [27] have 
shown the equivalence of the computation of the equilibrium dynamic structure factor of 
a Maxwell fluid both using fluctuating hydrodynamics and Mountain's approach. It is 
a major aim of the present paper to profit from the results of these previous works and 
to generalize them to the case where the fluid is subjected to an external temperature 
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gradient. The idea is to assess to what extent non-Newtonian effects, which have been 
considered negligible in actual experiments, do affect the light-scattering spectrum. 

The paper is organized as follows. In Sect. II we compute the equilibrium dynamic 
structure factor of a Maxwell fluid in which thermal conduction is governed by Fourier's 
law. This is done using fluctuating hydrodynamics. Sect. Ill deals with the same fluid 
which has now been brought to a steady state by the action of an externally imposed 
temperature gradient. Here we compute the light-scattering spectrum with the aid of the 
extension of fluctuating hydrodynamics to nonequilibrium steady states. The paper is 
closed in Sect. IV with a discussion and some concluding remarks. 

2 DYNAMIC STRUCTURE FACTOR FOR A MAXWEL 
FLUID IN EQUILIBRIUM. 

We here first analyze the case of a Maxwell fluid in equilibrium (i.e. a fluid in which 
instead of the Navier- Newton law, both elastic and viscous effects are incorporated into the 
constitutive equation for the stress tensor, c.f. Eq. (0) below) and we further assume that 
heat conduction is governed by the classical Fourier law. For this case, the corresponding 
full hydrodynamic equations (written in terms of the density, velocity and entropy density 
of the fluid denoted by p, v and s, respectively) are the continuity equation, the momentum 
balance equation and the energy balance equation, namely 

dp 

l£ + vVp = -pV-v, (!) 
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p— + pv-Vv = -Vp- V- r, (2) 

and 

ds <-» 
pT— + pTv-Vs = -V ■ q- r: Vv, (3) 
at 

which must be complemented with the constitutive equations for the heat flux vector q 
and the stress tensor r . These equations in our case read 



q = -AVT (4) 

and 

-r r ^- - r r v-V r = r +2t] Vv +rj v (V • v) 1 . (5) 
at 

Here, p is the hydrostatic pressure, T is the temperature, rj is the shear viscosity, 7] v is 
the bulk viscosity, C p is the constant density heat capacity, A is the thermal conductivity, 
1 is the unit tensor, Vv= ^ Vv + (Vv) — ^ V • v 1 stands for the symmetric and 
traceless part of the velocity gradient and r r is the Maxwell relaxation time. In Eq. 
(fj) one can eliminate the hydrostatic pressure in favor of p and s using the fact that 
Vp = — c 2 Vp + c^rVs (where c is the adiabatic speed of sound in the fluid, a is the 
thermal expansion coefficient and kj< is the isothermal compressibility) but the form we 



have chosen will prove convenient later on. Note also that on the left hand side of Eq. 



we have considered the usual material time derivative rather than the corrotational time 
derivative. Although the latter choice would be more satisfactory from the rheological 
point of view, Eq. (|5]) represents the simplest form of the Maxwell model. In any case, 
since later on we will retain only the linear terms, the use of the material time derivative 
will not imply loss of generality. 

We consider for simplicity the case where a = (which implies that C p = C p , where 
C p is the constant pressure heat capacity) and take r], r\ v , c, A, C p and r r to be constant. 
Let a subscript zero refer to the equilibrium value and a subscript one denote the (small) 
deviation from equilibrium, so that p = p + pi, v = Vi, s = Sq + S\, T = T + Tiand 
r = r 1 . Keeping only linear contributions in the deviations, one gets from Eqs. ([I]), and 



d ) to (1) 



dpi 



PoV • vi, 



(6) 



dt 



~dt 



1 



-c 2 Vpi - V- Ti, 



(7) 



PoT 
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and 
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-T r ^- = r 1 +2 V V° Vl +T] V (V • Vi) T . (9) 
at 

Within fluctuating hydrodynamics ||, stochastic contributions a and g are added to 
the momentum and heat fluxes, respectively. At this stage it is convenient to work in 
Fourier space. Denoting by an upper tilde the space and time Fourier transform of a 
given quantity, the linearized continuity, momentum and energy equations in (k, u) - 
space (with k the wavevector and u the frequency) are written as 



-iujpi = -ip it-vi, (10) 



~ 2i ~ (1 + iOJT r ) ( A \ , 2~ 

-zcjpovi = -ic kpi - — — —— [-T] + r] v )k vi - ik- a, (11^ 

1 + uj 2 t? \ 3 



and 



— iup T si = —\k 2 Ti — ik-g. (12) 

Note that in writing Eq. (|ll|) we have made use of the (k, uS) - space version of Eq. @. 
namely 



In order to close the fluctuating hydrodynamics approach, we now require the corre- 
sponding fluctuation-dissipation relations. Since in the (r, u)- space it follows from Eq. 
d) that 



1 + iuJT r ) 



2r] Vvi (r, (jj) + r) v (V • Vi (r, u)) 1 



1 + cu 2 r r 2 

2 V d (lu) V° Vl (r, u) + r)°(u) (V • vi (r, u)) T : 



(14) 



i.e. in the frequency domain the Maxwell model implies a linear relationship between 
the stress tensor and the thermodynamic forces but with dispersive transport coefficients, 
r] D (u) and r]^{ui) being the dispersive (complex and frequency-dependent) transport co- 
efficients, then following Landau and Lifshitz || the fluctuation-dissipation relations read 



(k,u;) #*(k', u')) = 2k B Tl (2tt) 4 \5ij5 3 (k - k') 5 (u - u') 



(15) 



and 



2k B T 
1 + ufir? 



(27T) 



2r]Aiji m + [ rj v - -r] ) <%<^ m 



5 3 (k - k') 5 (w - a/) 
(16) 



(ff^ (k,w)^*(k'V)) = (o« (k,uj)) = (gi(k,uj)) = 



(17) 
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the angular brackets denoting an (equilibrium) ensemble average over the initial states of 
the system, the asterisk meaning complex conjugation and Aiji m = | (Su5j m + 5i m 5ji), 6 rs 
being the Kronecker delta. 

Note that the energy equation within the present approximation is decoupled from 
both the continuity equation and the momentum balance equation. Since the dynamic 
structure factor S(k,uj) is related to the density-density correlation function through 



(pi (M pi (k>')> = (2tt) 4 5 3 (k - k') 5 (w - <J) S% (k,w) , (18) 

the subscript B having been introduced to stress the fact that, within the present ap- 
proximation, it will give the Brillouin spectrum, and since from Eqs. (fLOD , (pi]) and flT6|) 
it follows that 



(pi (k,w) p[ (kV)> = (2tt) 4 5 3 (k - k') 5 (w - J 



x- 



2A; s T (l + aA- 2 ) (fry + rj v ) k 4 



(a; 2 - c 2 A; 2 ) (1 + cjV 2 ) - ^ (§77 + 77J fc 2 J + g (| 77 + r) v ) k 4 
one readily identifies 



(19) 



(k,w) 



2fc B T (l + cu 2 r 2 )(|r ? + ^) A; 4 



- c 2 F) (1 + cu 2 r 2 ) - u 2 ^ (§77 + Tfc) fc 



(20) 
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The spectrum given in Eq. (|2"0| ) consists of two (symmetric) Brillouin peaks slightly 
displaced from u± = ±ck and a peak centered at uj c = 0, known as the Mountain peak, 



whose origin can be ascribed to the viscoelastic character of the ffuid[28 



Although in the approximation we have taken of a = 0, there is no Rayleigh line, if 
this restriction were avoided, the Rayleigh line, also centered around u c = 0, would follow 
from the entropy-entropy correlation function. In any case, since to a first approximation 
the Rayleigh line is determined by the coupling of the temperature fluctuations to the 
transverse velocity fluctuations [2(J, the missing terms (proportional to a) in Eqs. (|TT|) and 
(|T2| ) will play no role in the final result. So we compute the entropy-entropy correlation 
function by manipulating Eqs. (|T2|) and ([15]). The result is 



(si (k, u) s 1 (k, u)) = (k,u) = — , (21) 

Pq ~r LJrpK, ) 



where Dt = is the thermal diffusivity and the subscript R serves to make reference 



poc 



to the connection with the Rayleigh line. 

3 FLUCTUATING HYDRODYNAMICS AND THE 

LIGHT-SCATTERING SPECTRUM OF A MAXWELL 
FLUID UNDER AN EXTERNAL TEMPERATURE 
GRADIENT. 

We now examine the same Maxwell fluid in which heat conduction is governed by Fourier 's 
law but in a layer confined between two parallel plates at z — +|f and z = — and of 
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infinite extension in the xy— plane. In the fluid layer a steady temperature gradient is 
maintained by fixing the temperatures of the upper and lower planes at different values, 
namely T% at z = — -| and T 2 at z = +-|. We assume that neither the gravitational force 
nor any other external forces act on the system. A time-independent solution to Eqs. ([]]) 
to (Q) leading to a stationary state (denoted by a subscript s) for the system satisfies the 
relations 



v s = (22) 



V Ps = (23) 



V • q s = 0. (24) 
Due to the symmetry of the problem, Eq. (f24"|) together with Eq.(HI) leads to 



For simplicity, let us look at the case where a s = and take rj s , rj vs , c s , X s , C ps and 
r rs to be constant. Since X s has been assumed to be a constant, it follows from Eq. 
and the boundary conditions that 

T s = T - Az 
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where 



rp _|_ rp 

To = = T S (Q) 



and 

Lj 

Further, the condition implied by Eq. (|23|) leads to Vp s = or p s = p 0; while Eq. 
( |22| ) implies that t s =0. Let a subscript Is now denote the (small) deviation from the 
steady state, so that p = p s + pu, v = v s + v ls = v ls , s = s s + si s , T = T s + T ls , 
q = q s + qi s = — A S VT S + qi s and t = t s + r ls =r ls . Keeping only linear contributions 
in the deviations, one gets from Eqs. (|J) to (H) 



dpu 
dt 



-p s V • v ls , (26) 



Ps^f= -c^Pis - V- r ls , (27) 



at 



q ls = -A s VT ls , (29) 
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and 



-t ts dt s = r u +2r] s Vvi s +r] vs (V • v is ) 1 . (30) 



The balance equations in the (r,u)— domain read 



-iupu = -p s V • vis, (31) 



icvp s v u = -c s Vp ls - V- ri,, (32) 



and 



-iup T s ls + Ap C ps Vi s • z = -V • q ls . (33) 

where z is a unit vector in the z direction. We now add fluctuating contributions to the 
heat flux and stress tensor denoted by g s and a s , respectively, take the spatial Fourier 
transform and make use of the stationary state analog of Eq. (ITJ) to arrive at 



-iup u = -ip s k ■ vi a , (34) 



■ ~ • 2 , ~ ( 4 / 3 ^ + Vvs) (1 + wr rs /c 2 _ <-> 
-2u;p s v ls = -zc s kp ls t— — - «k • <r s , (35) 

U + ^ C rJ 

13 



and 



-iup T s ls = -p T D Ts k 2 s ls -ik-g s - Ap C ps v ls ■ z. 



(36) 



In order to proceed any further, we require the fluctuation-dissipation relations obeyed 
by g s and a s . Under the hypothesis mentioned in the Introduction, we follow previous 
authors in assuming that the fluctuation-dissipation relations have the same form as that 
for thermal equilibrium, but with the equilibrium temperature replaced by the stationary 
temperature and the average taken over an stationary ensemble. Thus, restricting to 
small gradients (linear order in A) in our case one gets 

(g is (k,cu) g* s (k', u')) = 2k B Tl (2tt) 4 X^S (u - J) 



x 



5 3 (k _ k ' ) + l AL (s 3 (k - k' + Z) - 5 3 (k - k' - - 



(37) 



M r s im(V, "0> = (2vr) 4 



2r} s A ijlm + ( r} vs - -r) s ) bijb im 



S(uj- J) 



while 



x 



,3 (k _ k0 + ^3A k , + ! 



<5 3 k 



k>-* 
L 



(38) 



(a sij (k,u) g* sl (k',u')) = (a sij (k,uj)) = (g sl (k,uj)) = 0. 



(39) 
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Tremblay, Arai and Siggia[TT] have shown that, in a light scattering experiment in 
which a fluid in a stationary state due to an external temperature gradient is involved, 



S B (k,u) = (p 1 (k,u)pl(k,u)) + (p 1 (k+z/2L,u)p* 1 (k-z/2L,u;)) 
+ (pt Qs.-z/2L,w)ft(\sL+z/2L,w)) 

is the correlation function of interest for the Brillouin spectrum. By a simple manipulation 
of Eqs. (|34"D , ( |35"D and (|38|) it is easy to arrive at 



S B (k, u) = S e B q (k, u 



1 - 



2c/ As (u) A (k - z) 



T [(u 2 - C 2 (w) fc 2 ) 2 + (w) fc 2 ) 2 ] 



fc B T p A s (g) fc 2 
C? (w) 2 



l-e(k,o;) 



+ 



l+e(k,o;) 



(w - C s (w) A;) z + i (As (w) fc 3 )* (w + C s (w) fc) z + i (A s (w) fc 2 )' 

(40) 



where S B 9 (k, to) is given by Eq. fl20"D and we have identified p s = po,?7s = ?7, ?7us = T) v 
c a = c, and r rs = r r and introduced the following quantities 



and 



A, (u) 



A>(1 + <A?)' 



C s (w) = [c 2 + LU 2 T r D ls (w)] 



1/2 



e(k,w) 



C s (w) (k • z) 
T As H k 2 



2 (loD 1s (u) k 2 f 



(lo 2 -C 2 (oo) k 2 ) 2 + (uD ls (oo)k 2 ) 2 



(41) 
(42) 



(43) 
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It should be pointed out that in order to perform the expansion leading to the approxi- 
mation appearing in the second line of Eq. fllPf) , apart from the weak gradient assumption 
(i. e. ^4^- *C 1) it has also been assumed that 

Now we turn to the Rayleigh line. Again the restriction a s = implies that the 
spectrum in Eq. (pE0|) contains no Rayleigh line. If such restriction is removed, in order 
to compute the Rayleigh line one again requires the entropy-entropy correlation function. 
By the same line of reasoning used in connection with Eq. ( ]2lD , the missing terms 
(proportional to a s ) in Eqs. ([35]) and (|36|) will play no role in the final result. So we 
proceed with the calculation. From Eq. fl36|) and making the identification Dt s = it 
follows that 



c s M(k-g) 



< 1. 



~ /, \ ~ ik ■ gj (PqTq) - ACpV ls ■ z/Tp 

s ls k, u) = ; — 2 • 44 

— iuj + DtK 

We only need the transverse velocity fluctuations v ts which are obtained from v ls as 



v ts = v ls • ( 1 - kkV z. (45) 



It is clear from this last equation that if the scattering vector k is taken parallel to 
the temperature gradient, then v ts = 0, so there will be no effect on the entropy-entropy 
correlation function. Therefore we choose k _L z. From Eqs. fl35|) and (|I"4]) it follows that 



zk • a s ■ ( 1 — kk) • z 



~ iUJ P0 + (lW) 
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Using this expression in Eq. yields 



Si, (k, u) = —=-. —= — + ) r . (47) 

~'~ + Po(lW^) 



. iACpk • cr s • ^ 1 - kkj • z 

A,T (-ia; + D T k 2 ) pQ T Q {-iu + D T £; 2 ) (-i 

Since any asymmetries occur in the z direction, in the scattering plane (which is 
orthogonal to the temperature gradient) spatial and temporal invariance holds. This in 
turn implies that from the fluctuation-dissipation relations one finds (for k, k' _L z) 

(k ■ g s (k,w) k' • g:(k', u')) = 2k B T 2 (2vr) 4 Po C p D T k 2 6 3 (k - k') 6 (u - cu') , (48) 

and 



k • a s ■ (T - kk) • z • k' ■ a] ■ (T - k'k') ■ za^ = (2tt) 4 V k 2 5 3 (k - k') 5 (u - to') 

(49) 

Use of Eqs. © to © leads finally to 

(Si, (k, w) s* ls (k, w)) = (k, u) 



a 



1 + 



Pp (l+U> 2 T?) 



A 2 



T D T 



u 2 1 + 



+ 



r; 2 fc 4 



(k 



(50) 



with S 1 ^ (k, uj) given by Eq. (pT|). Eqs. (f40"|) and ( p0|) are the main results of this paper. 
A brief assessment of their potential use will be carried out in the following Section. 
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4 DISCUSSION. 

As a first observation, we note that if the limit r r — * is taken, then Eqs. (p0[), (pip , 
PD| ) and fl5*0"| ) reduce to the well known results for the Newtonian fluid (c.f. for instance 
Refs. |5|, [|ll|] and f2"0fl), as they should. Also, the equilibrium light-scattering spectrum 



in this paper, Eq. (|20|), arises as a particular case of the one computed in Ref. 
provided in the latter the relaxation times of the trace and traceless symmetric part of 
the stress tensor are both set equal to r r , the internal degrees of freedom are eliminated 
and the approximation a = is made. In order to get a picture of the importance of the 
viscoelastic effects in equilibrium, in Fig. 1 we display the equilibrium spectrum (with 
a = 0) of both the Maxwell fluid and the Newtonian fluid. The thermophysical properties, 
given in Table 1, are those of toluene |29[| , the temperature has been set to T = 277 K, 
we have taken as wavenumber k = 85200 vrC x (which is the value corresponding to a 

o 

scattering angle of 10 mrad and a wavelength of 4880 A) and we have included a very 
crude estimate of the relaxation time r r ~ 10~ 5 sec of this liquid obtained from Grad's 
kinetic theory expression. 

Concerning the nonequilibrium spectrum again with a = 0, the presence of the external 
temperature gradient has a similar effect to the one observed in the Newtonian fluids 
for the Brillouin peaks, as illustrated in Fig. 2. Here we have again used the data 
corresponding to toluene and taken the gradient to point upwards with A = 100 K/m, 
L = 0.01 m and T = 277 K, these latter numbers fulfilling the restrictions involved 
in our derivation. A noteworthy feature is that the asymmetry in the Brillouin peaks, 
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related to the fact that there are more sound waves absorbed with wavevector — k than 
those emitted with wavevector k due to the temperature gradient and its corresponding 
heat flux is much more pronounced and the peaks become narrower for the viscoelastic 
fluid. On the other hand, within the various approximations and simplifying asumptions 
that we made, the Mountain peak is not affected by the temperature gradient. This, 
which is somewhat surprising and most likely tied to the assumption a = 0, is a trivial 
consequence of the fact that the nonequilibrium correction vanishes for u c = 0. 

In conclusion, we want to point out that we have derived (to our knowledge for the 
first time) the light-scattering spectrum for a Maxwell fluid in a steady state under the 
action of an external temperature gradient. Although many (reasonable) approximations, 
some of which might be easily disposed of, were made along the way in order to get 
analytical results, and in spite of the fact that the Maxwell fluid is a very simple model of 
a viscoelastic fluid, we are persuaded that the present calculation captures the essential 
(qualitative) features to be observed in light-scattering experiments in connection with 
the elasticity of fluids. Whether this expectation is real awaits the necessary confrontation 
with experiment. 
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Figure Captions. 

Fig. 1 Equilibrium light-scattering spectrum as a function of frequency u for a fluid 
whose thermophysical data correspond to those of toluene at T = 277 K. The scattering 
wave-number is k = 85200 mr l . The dashed line represents the results for the Newtonian 
fluid and the continuous line those of the Maxwell fluid. 

Fig. 2 Light-scattering spectrum as a function of frequency uj for a fluid whose thermo- 
physical data correspond to those of toluene at T = 277 K under the action of an external 
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temperature gradient (A = 100 K/m and L = 0.01 K). The scattering wave-number is 
k = 85200 m~ l . The code is the same as that of Fig. 1. 

Table 1 

rj = 0.000553 Pa. sec 
r] v = 0.00747 Pa. sec 
p = 861.5 Kg/m 3 
c = 1285 mj sec 

Table Caption. 

Table 1. Thermophysical data of toluene as given in Ref. 
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